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Abstract: We consider the standard first passage percolation model in the rescaled graph l//n 

\ for d > 2, and a domain Q. of boundary T in W^. Let and be two disjoint open subsets 

^ ■ of r, representing the parts of T through which some water can enter and escape from Q.. We 

10 ■ investigate the asymptotic behaviour of the flow (j)n through a discrete version 0„ of O between 

\ the corresponding discrete sets T\ and T^. We prove that under some conditions on the regularity 

I of the domain and on the law of the capacity of the edges, converges almost surely towards 

' a constant 4>^^ which is the solution of a continuous non-random min-cut problem. Moreover, we 

Q , give a necessary and sufficient condition on the law of the capacity of the edges to ensure that </>f7 > 0. 
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1 First definitions and main result 

We use many notations introduced in [18] and [H]. Let d > 2. We consider the graph (Z^,E^) 
having for vertices = Z'^/n and for edges E^, the set of pairs of nearest neighbours for the 
standard norm. With each edge e in we associate a random variable t{e) with values in 
M+. We suppose that the family (t(e),e G E^) is independent and identically distributed, with a 
common law A: this is the standard model of first passage percolation on the graph (Z5(,E5(). We 
interpret t{e) as the capacity of the edge e; it means that t{e) is the maximal amount of fluid that 
can go through the edge e per unit of time. 

We consider an open bounded connected subset 17 of such that the boundary F = dQ of 17 
is piecewise of class (in particular F has finite area: 7{'^~^{T) < 00). It means that F is included 
in the union of a finite number of hypersurfaces of class C^, i.e., in the union of a finite number of 

submanifolds of M'' of codimension 1. Let F^, F^ be two disjoint subsets of F that are open in F 
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We want to define the maximal flow from to through Q for the capacities (t(e), e € EJ^). We 
consider a discrete version {0,n,^n,^n,^n) of (^)r,r^,r^) defined by: 

= {xe Zf^\doo{x,n) < 1/n}, 
r„ = {x enn\3y ^^n, {x, y) e E^} , 

= {x ern\doo{x,r) <l/n, doo{x,r^-')>l/n} for i = l,2, 

where doo is the L°°-distance, the notation {x,y) corresponds to the edge of endpoints x and y (see 
figure [I|) . 




Figure 1: Domain 0,. 

We shall study the maximal flow from to in O^. Let us define properly the maximal 
flow 4>{Fi — > F2 in C) from Fi to F2 in C, for C C M'^ (or by commodity the corresponding graph 
C n Z'^/n). We will say that an edge e = {x, y) belongs to a subset A of M"^, which we denote by 
e & A, if the interior of the segment joining x to y is included in A. We define E^ as the set of all the 
oriented edges, i.e., an element e in Kf^ is an ordered pair of vertices which are nearest neighbours. 
We denote an element e G E^ by {{x, y)), where x, ?/ G are the endpoints of e and the edge is 
oriented from x towards y. We consider the set S of all pairs of functions {g, o), with g : E^ — > R"*" 
and o : E^ — > EJ^ such that o{{x, y)) e {((x, y)), ((y, x))}, satisfying: 

• for each edge e in C we have 

< g{e) < t{e) , 

• for each vertex f in C \ (Fi U F2) we have 

eeC : oie)=((v,-)) eeC : o(e)=((-,v)) 

where the notation o(e) = {{v,.)) (respectively o(e) = {{.,v))) means that there exists y G such 
that e = {v,y) and o(e) = {{v,y)) (respectively o(e) = {{y,v))). A couple (ff, o) € 5 is a possible 
stream in C from Fi to F2: g{e) is the amount of fluid that goes through the edge e, and o(e) gives 
the direction in which the fluid goes through e. The two conditions on (g, o) express only the fact 
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that the amount of fluid that can go through an edge is bounded by its capacity, and that there is 
no loss of fluid in the graph. With each possible stream we associate the corresponding flow 

flow(g,o) = ^ 9((^^,'f^»lo((«,«»=((«,«» - S'((^^,^^>)lo((n,i,»={(«,n» ■ 

This is the amount of fluid that crosses C from Fi to F2 if the fluid respects the stream {g,o). The 
maximal flow through C from Fi to F2 is the supremum of this quantity over all possible choices of 
streams 

(P{Fi ^ F2 in C) = sup{flow(g, o) I (g, o) e S} . 

We denote by 

<An = 4>{ri ^ r2 in nn) 

the maximal flow from to in Qn- We will investigate the asymptotic behaviour of (j)n/n'^~^ 
when n goes to infinity. More precisely, we will show that {4>n/n'^~^)n>i converges towards a 
constant (f)n (depending on $7, F^, F^, A and d) when n goes to infinity, and that this constant is 
strictly positive if and only if A(0) < 1 — pc{d), where Pc{d) is the critical parameter for the bond 
percolation on Z*^. The description of cf)Q will be given in section [2l Here we state the precise 
theorem: 

Theorem 1. We suppose that Q is a Lipschitz domain and that F is included in the union of a 
finite number of oriented hypersurfaces Si, ...,Sr of class which are transverse to each other. We 
also suppose that F^ and F^ are open in T, that their relative boundaries drT^ and dpT'^ in F have 
null 7i'^~^ measure, and that (i(F^,F^) > 0. We suppose that the law A of the capacity of an edge 
admits an exponential moment: 

39 >0 / e^'^'dA^x) < +00 . 
Jr+ 

Then there exists a finite constant (pQ > such that 

1™ = a.s. 

n—^oo n 

Moreover, this equivalence holds: 

0n > ^ A(0) < 1 - p,{d) . 

Remark 1. In the two companion papers [7] and [8], we prove in fact that the lower large deviations 
of 4>n/n'^~^ below (pQ are of surface order, and that the upper large deviations of (j)n/n'^~^ above (po, 
are of volume order (see section [3^ where these results are presented). 

2 Computation of 4)^ 
2.1 Geometric notations 

We start with some geometric definitions. For a subset X of M"^, we denote by 1-L^{X) the s- 
dimensional Hausdorff measure of X (we will use s = d — 1 and s = d — 2). The r-neighbourhood 
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Vi{X,r) of X for the distance d,, that can be the Euchdean distance if z = 2 or the L°°-distance if 
i = oo, is defined by 

ViiX,r) = {yeM.''\d,{y,X)<r}. 

If X is a subset of included in an hyperplane of M'* and of codimension 1 (for example a non 
degenerate hyperrectangle), we denote by hyp{X) the hyperplane spanned by X, and we denote by 
cyl(X, h) the cylinder of basis X and of height 2h defined by 

cyl(X, h) = {x + tv\x e X ,t e [-h, h]} , 

where v is one of the two unit vectors orthogonal to hyp(X) (see figure [2]). For x € M'^, r > 




Figure 2: Cylinder cyl{X,h). 

and a unit vector u, we denote by B{x,r) the closed ball centered at x of radius r, by disc(x,r, u) 
the closed disc centered at x of radius r and normal vector v, and by hyp{x,v) the hyperplane 
containing x and orthogonal to v. We denote by the volume of a unit ball in M"^, and a^-i the 
7i'^~^ measure of a unit disc. 

2.2 Flow in a cylinder 

Here are some particular definitions of flows through a box. Let A be a non degenerate hyperrect- 
angle, i.e., a box of dimension d — 1 in W^. All hyp err ect angles will be supposed to be closed in W^. 
We denote by v one of the two unit vectors orthogonal to hyp(74). For h a positive real number, we 
consider the cylinder cyl{A,h). The set cyl{A,h) \ hyp(^) has two connected components, which 
we denote by Ci{A, h) and C2{A, h). For i = 1, 2, let A^ be the set of the points in Ci{A, h) n 
which have a nearest neighbour in Zj^ \ cy\{A, h): 

A\ = {x G Ci(^, /i) n I 3y G \ cyl(A, h) , (x, y) G EJ^} . 

Let T(A, h) (respectively i?(A, /i)) be the top (respectively the bottom) of cyl(yl, /i), i.e., 

T[A, h) = {x £ cyl{A, h)\3y ^ cy\{A, h) , (x, y) G and (x, y) intersects A + hv} 

and 

B{A, h) = {x e cyl{A, h)\3y ^ cyl{A, h) , {x, y) G E^ and (x, y) intersects A - hv] . 
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2.3 Max-Bow min-cut theorem 



For a given realisation (t(e),e G E,„) we define the variable T{A,h) = T{cyl{A, h),v) by 

TiA,h) = TicyliA,h),v) = </.(4 ^4 in cyliA,h)), 

and the variable (j){A,h) = (/)(cyl(A, /i), f ) by 

cP{A,h) = cl)icy\{A,h),v) = <j)iBiA,h)^T{A,h) in cyliA,h)), 

where (/>(-Fi F2 in C) is the maximal flow from Fi to F2 in C, for C C M'^ (or by commodity the 
corresponding graph C PlZ'^/n) defined previously. The dependence in n is implicit here, in fact we 
can also write r„(A, /i) and if we want to emphasize this dependence on the mesh of the 

graph. 



2.3 Max-flow min-cut theorem 

The maximal flow (piFi — > F2 in C) can be expressed differently thanks to the max-flow min-cut 
theorem (see [5]). We need some definitions to state this result. A path on the graph from vq 
to Vjn is a sequence {vo,ei,vi, ...,em,Vm) of vertices vq, ...,Vm alternating with edges ei, such 
that Vi-i and Vi are neighbours in the graph, joined by the edge e,, for i in {1, ...,m}. A set E of 
edges in C is said to cut Fi from F2 in C if there is no path from Fi to F2 in C \ E. We call E an 
(Fi, F2)-cut if E cuts Fi from F2 in C and if no proper subset of E does. With each set E of edges 
we associate its capacity which is the variable 

V{E) = Y,t{e). 

The max-flow min-cut theorem states that 

(p{Fi F2 in C) = min{ V{E) | F is a (Fi, F2)-cut } . 
In fact, as we will see in section [231 (pQ is a continuous equivalent of the discrete min-cut. 



2.4 Definition of u 

The asymptotic behaviour of the rescaled expectation of Tn{A, h) for large n is well known, thanks 
to the almost subadditivity of this variable. We recall the following result: 

Theorem 2. We suppose that 

I X dA{x) < 00 . 

-'[0,4- oo[ 

Then for each unit vector v there exists a constant u(d,A,v) = I'iv) (the dependence on d and A is 
implicit) such that for every non degenerate hyperrectangle A orthogonal to v and for every strictly 
positive constant h, we have 

r nrn{A,h)] 



+00 n' 



5 



2.5 Continuous niin-cut 



2 COMPUTATION OF (pn 



For a proof of this proposition, see [25]. We emphasize the fact that the hmit depends on the 
direction of v, but not on h nor on the hyp err ect angle A itself. 

We recall some geometric properties of the map u : v 5^^"^ i— > i^iv), under the only condition 
on A that E(t(e)) < oo. They have been stated in the section 4.4 of [25]. There exists a unit vector 
vq such that I'ivo) = if and only if for all unit vector v, u{v) = 0, and it happens if and only 
if A({0}) > 1 — Pcid). This property has been proved by Zhang in [27]. Moreover, v satisfies the 
weak triangle inequality, i.e., if (ABC) is a non degenerate triangle in M*^ and va, vb and vc are 
the exterior normal unit vectors to the sides [BC], [AC], i-^B] in the plane spanned by A, B, C, 
then 

■H\[AB])v{vc) < nH[AC])iy{vB) +nH[BC]MvA) ■ 
This implies that the homogeneous extension uq of to M'^, defined by fo(0) = and for all w in 

uoiw) = \w\2iy{w/\w\2) , 

is a convex function; in particular, since is finite, it is continuous on W^. We denote by fmin 
(respectively Vmax) the infimum (respectively supremum) of v on S'^~^. 

2.5 Continuous min-cut 

We give here a definition of 4>q and of another constant cpn in terms of the map u. For a subset F 
of M"', we define the perimeter of -F in by 

V{F,n) = sup|y^div/(x)d£'^(x), /eC-(f^,i?(0,l))| , 

where C^i^, B{0, 1)) is the set of the functions of class C°° from to B{0, 1), the ball centered at 
and of radius 1 in W^, having a compact support included in and div is the usual divergence 
operator. The perimeter V{F) of F is defined as V{F,M.'^). We denote by dF the boundary of F, 
and by d*F the reduced boundary of F. At any point x of d*F, the set F admits a unit exterior 
normal vector vpix) at x in a measure theoretic sense (for definitions see for example [9] section 
13). For all F C M"' of finite perimeter in we define 

In{F) = [ y{vF{x))dn''-\x)+ f iy{v(^Fr.^){x))dn''-\x) 
Jd*Fr\n Jr2n9*(Fnn) 

+ / iy{vn{x))dn'^-\x) . 

^rin9*(n\F) 

If 'P(F, Q) = +00, we define In{F) = +oo. Finally, we define 

(l)n = mf{IniF) \F cM'^} = inf{Jn(F) {FcQ}. 
In the case where dF is C^, 1q,{F) has the simpler following expression: 

In{F) = [ u{vF{x))dn''-\x)+ [ iy{v(^Fr.^){x))dn''-\x) 
JdFnn Jr^ndiFnn) 

+ / i^Mx))dn'^-\x) . 
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2.5 Continuous niin-cut 




{dF nn)u (r^ n d{F nn)) u (r^ n d{Q \ F)) 
Figure 3: The set {dF nn)U {T^ n d{F n Q)) U (F^ n 9(0 \ F)). 



The locahzation of the set along which the previous integrals are done is illustrated in figure [HI 

When a hypersurface S is piecewise of class C^, we say that S is transverse to T if for all x S 5nr, 
the normal unit vectors to S and T at x are not collinear; if the normal vector to S (respectively 
to r) at X is not well defined, this property must be satisfied by all the vectors which are limits of 
normal unit vectors to S (respectively T) at y £ S (respectively y £T) when we send y to x - there 
is at most a finite number of such limits. We say that a subset P ofW^ is polyhedral if its boundary 
dP is included in the union of a finite number of hyperplanes. For each point x of such a set P 
which is on the interior of one face of dP, we denote by vp{x) the exterior unit vector orthogonal 

o 

to P at x. For A C M"', we denote by A the interior of A. We define 4>q by 



inillniP) 



P C M"^ , Fi C P , F2 c M'' \ P 
P is polyhedral , dP is transverse to F 



Notice that if P is a set such that 



Fi C P and F2 c R'^ \ P , 



then 



MP) 



dPnn 



v{vp{x))<m'^-^{x). 



See figure m to have an example of such a polyhedral set P. 

The definitions of the constants 4)0, and cf)^ are not very intuitive. We propose to define the 
notion of a continuous cutset to have a better understanding of these constants. We say that 5 C M'^ 
cuts F^ from F^ in O if every continuous path from F^ to F^ in O intersects S. In fact, if P is a 
polyhedral set of such that 



Fi C P and F2 c M'^ \ P , 

then dP H is a continuous cutset from F^ to F^ in 0. Since v{v) is the average amount of fluid 
that can cross a hypersurface of area one in the direction v per unit of time, it can be interpreted 
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dP 



/ ^ on 




Figure 4: A polyhedral set P as in the definition of 



as the capacity of a unitary hypersurface. Thus 1^{P) can be interpreted as the capacity of the 
continuous cutset S-PnO. The constant (f)^ is the solution of a min cut problem, because it is equal 
to the infimum of the capacity of a continuous cutset that satisfies some specific properties. We 
can define two other constants, that are solutions of possibly more intuitive min cuts problems. If 
5 is a hypersurface which is piecewise of class C^, we denote by vs{x) one of the two normal unit 
vectors to 5 at x for every point x at which S is regular. The measure of the points at which 

S is not regular is null. We define 



inf<^ / _v{vs{x))(m'^-^{x) 
Jsnn 



S hypersurface piecewise of class 
S cuts from in Q 



and 



<Pn = inf <^ / _uivsix))dn''^'ix) 
Jsnn 



S polyhedral hypersurface 
S cuts F^ from F^ in Q 



We remark that by definition, 

(t>n < 4>Q < (t>n- 

We claim that (pn ^ 4>n- Let 5 be a hypersurface which is piecewise of class C^, which cuts F^ from 
F^ in 0, and such that 

/ _u{vs{x))dn''-\x) <^+r] 
Jsnn 

for some positive rj. Let F be the set of the points of O \ 5 that can be joined to a point of F^ by 
a continuous path. Then 

{dFnn)u{r^nd{n\F))u{r'^nd{Fnn)) c snTi. 

Thus F is of finite perimeter in O, and 1n{F) satisfies 

In{F) < I _v{vs{x))dH''~\x) <^ + rj. 
Jsnn 
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Thus we have proved that ^ 

3 State of the art 

3.1 Existing laws of large numbers 

Only in this section, we consider the standard first passage percolation model on the graph (Z'^, E*^) 
instead of the rescaled graph (Z^,E5^). We present here some laws of large numbers that have been 
proved about maximal flows. 

Using a subadditive argument and concentration inequalities, Rossignol and Theret have proved 
in [25] that T{nA, h{n)) satisfies a law of large numbers: 

Theorem 3 (Rossignol and Theret). We suppose that 

I X dK{x) < oo . 

i[0,oo[ 

For every unit vector v, for every non degenerate hyperrectangle A orthogonal to v, for every height 
function /i : N ^ satisfying lim„,^oo h{n) = +oo, we have 

T{nA,h{n)) ^ 

Moreover, if the origin of the graph belongs to A, or if 

I x^+d^ dK{x) < oo , 

i[0,oo[ 

then 



T{nA,h{n)) 



n— >oo 



Kesten, Zhang, Rossignol and Theret have studied the maximal flow between the top and the 
bottom of straight cylinders. Let us denote by D{k, m) the cylinder 



d-i 



D{k,m) = Y[[0,ki] X [0, 



i=l 



where k = (fci, ...,kd-i) G M'^ ^. We denote by ^(k, m) the maximal flow in i3)(k, m) from its top 
nf=i [0' ^i] ^ {"^} to its bottom Ilti [0' ^i] ^ {0}- Kesten proved in [19] the following result: 

Theorem 4 (Kesten). Let d = 3. We suppose that A(0) < po for some fixed po > 1/27, and that 

37 > / e^^ dA{x) < 00 . 

If m = m(k) goes to infinity with ki > ^2 in such a way that 

36 >0 lim A;"^+'^logm(k) = 0, 
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then 

(/)(k, TTl) 1 

lim — — = z^((0,0, 1)) a.s. and in L . 

ki>k2^oo kik2 

Moreover, if A{0) > l—pc{d), where Pc{d) is the critical parameter for the standard bond percolation 
model on Z"^, and if 

I dk{x) < oo, 

-'[0,+oo[ 

there exists a constant C = C{F) < oo such that for all m = m(k) that goes to infinity with ki > /c2 
and satisfies 

. m(k) 
hmmf ------ > C, 

ki>k2~*oo kik2 

for all ki > k2 sufficiently large, we have 

(^(k, m) = a.s. 
Zhang improved this result in |28| where he proved the following theorem: 
Theorem 5 (Zhang). Let d>2. We suppose that 

37 > / e^^ dk{x) < 00 . 

-'[0,+oo[ 

Then for all m = m(k) that goes to infinity when all the ki, « = 1, d — 1 go to infinity in such a 
way that 

35g]0,1] logm(k) < max k]~\ 

j=l,...,d— 1 

we have 

lim — -7^, — = z^((0, 0, 1)) a.s. and in . 

fci,...,fcd-i-oo nf=i 

Moreover, this limit is positive if and only if A{{)) < 1 — Pc{d). 

To show this theorem, Zhang obtains first an important control on the number of edges in a 
minimal cutset. Finally, Rossignol and Theret improved Zhang's result in [25] in the particular case 
where the dimensions of the basis of the straight cylinder go to infinity all at the same speed. They 
obtain the following result: 

Theorem 6 (Rossignol and Theret). We suppose that 

X dA{x) < 00 . 



I 

J\o. 



'[0,oo[ 

For every straight hyperrectangle A = Jlti I^) ^i] ^ {0} a, > for all i, for every height function 
h -.N M"*" satisfying hm„_»oo h{n) = +00 and lim^^oo log h{n)/n'^~^ = 0, we have 

(j)(nA,h(n)) 1 
^™ o/^ 1/ = i^((0, •••,0, 1)) a.s. and m L . 

n^oo TL'^~'^[nA) 
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In dimension two, more results are known. We present here two of them. Rossignol and Theret 
have studied in [21] the maximal flow from the top to the bottom of a tilted cylinder in dimension 
two, and they have proved the following theorem (Corollary 2.10 in [24]): 

Theorem 7 (Rossignol and Theret). Let v be a unit vector, A a non degenerate line-segment orthog- 
onal io f , /i : N — > M'^ a height function satisfying lim„_^oo h{n) = +oo and lim„_>oo log h{n)/n = 0. 
We suppose that there exists a G [0,7r/2] such that 



2h(n) 
lim — j- — — 

n^oo n^(nA) 



tan a . 



Then, if 



we have 



[0,oo[ 



X dA{x) < oo , 



(f>{nA,h{n)) 



lim ^ , ^ , 



inf / — — \ v satisfies v ■ v' > cos a\ in 
yv ■ v' J 



Moreover, if the origin of the graph is the middle of A, or if 

dA{x) < oo , 



J [0,00 [ 



then we have 



(p{nA,h{n)) 
lim — — r- — — — 

n^oo n^{nA) 



inf 



V ■ V 



v' satisfies v ■ v' > cosaj' 



a.s. 



Garet studied in [15] the maximal flow cr^A) between a convex bounded set A and infinity in the 
case d = 2. By an extension of the max flow - min cut theorem to non finite graphs, Garet proves 
in [15] that this maximal flow is equal to the minimal capacity of a set of edges that cuts all paths 
from A to infinity. Let dA be the boundary of A, and d*A the set of the points x G dA at which 
A admits a unique exterior normal unit vector VAix) in a measure theoretic sense (see [9], section 
13, for a precise definition). If ^ is a convex set, the set d*A is also equal to the set of the points 
X £ dA at which A admits a unique exterior normal vector in the classical sense, and this vector is 
VAix). Garet proved the following theorem: 

Theorem 8 (Garet). Let d = 2. We suppose that A(0) < 1 - pc{2) = 1/2 and that 



37 > / e^^ dA{x) < 00 . 

J\0,-\-oo\ 



'[0,+oo 

Then for all convex bounded set A containing in its interior, we have 

a{nA) 



lim 



u{vA{x))dn\x) = I{A) > a.s. 



n— >oo n Jd*A 

Moreover, for all e > 0, there exist constants Ci, C2 > depending on e and A such that 

'a{nA) 



Vn > 



nI{A) 



\l-eA+e 



< Ci exp(— C2n) 
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Nevertheless, a law of large numbers for the maximal flow from the top to the bottom of a tilted 
cylinder for d > 3 was not proved yet. In fact, the lack of symmetry of the graph induced by the 
slope of the box is a major issue to extend the existing results concerning straight cylinders to tilted 
cylinders. The theorem of Garet was not extended to dimension d > 3 either. Theorem [1] applies 
to the maximal flow from the top to the bottom of a tilted cylinder. Thus it is a generalisation 
of the laws of large numbers of Kesten, Zhang, Rossignol and Theret for the variable 4> in straight 
cylinders, in the particular case where all the dimensions of the cylinder go to infinity at the same 
speed (or, equivalently, the cylinder is fixed and the mesh of the graph go to zero isotropically). 
Moreover, it gives a hint of what could be a generalisation of the result of Garet in higher dimension, 
all the more since the expression of the constant (pQ is a reminiscent of the value of the limit in 
Garet's Theorem: the capacity Iq of a continuous cutset is exactly the same as the one defined by 
Garet in [15] in dimension two, except that we consider a maximal flow through a bounded domain, 
so our capacity is adapted to deal with specific boundary conditions. 

Prom now on, we work in the rescaled graph (ZJ^, EJ^). 

3.2 Large deviations for 0„ 

We present here the two existing results concerning We consider an open bounded connected 
subset $7 of M*^, whose boundary F is piecewise of class C^, and two disjoint open subsets and 
F^ of F. The first result states that the lower large deviations below cpQ are of surface order, and is 
proved by the authors in [7]: 

Theorem 9. // the law A of the capacity of an edge admits an exponential moment: 

39 >0 / e^'''dA{x) < +oo , 
Jr+ 

and if A(0) < 1 — Pc{d), then for all A < 

hmsup --^ log P An''^^] < 0. 

n— too IT' 

The second result states that the upper large deviations of tpn above (pn are of volume order and 
is proved by the authors in [8]: 

Theorem 10. We suppose that d{T^,T'^) > 0. // the law A of the capacity of an edge admits an 
exponential moment: 

36 >0 [ e^^(iA(x) < +oo , 
Jr+ 

then for all \> (f)^, 

lim sup — log W[(j)n > \n'^~^] < . 

n— too 

By a simple Borel-Cantelli lemma, these results imply that if A admits an exponential moment 
and if d(F\F2) > 0, then 

(f)n < lim inf — < hm sup — < (f>n . 

n^oo n n^oo n 
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•n 



Notice here that Theorem [9] allows us to obtain the first inequality only under the additional 
hypothesis that A(0) < 1 — pc{d), however if A(0) > 1 — Pc{d) we know that ^{v) = for all u, so 
4)n = ^ and the first inequality remains valid. 

Thus, to prove Theorem [H it remains to prove that 4>n = 4>n, and to study the positivity of 
<pQ. The equality 4>q = is a consequence of a polyhedral approximation of sets having finite 
perimeter that will be done in section [H The positivity of (pQ is proved in section [H using tools of 
differential geometry like tubular neighbourhood of paths. These two results are proved by purely 
geometrical studies. Since the probabilistic part of the proof of Theorem [T] is contained in Theorems 
[QlandfTOl we propose a sketch of the proofs of these two theorems in sections [3.2. II and 13.2.21 to help 
the understanding of the law of large numbers proved in this paper. 

Before these two sketches of proofs, we would like to make two remarks. The first one is that the 
large deviations that are obtained in Theorem [9] and [10] are of the relevant order. Indeed, if all the 
edges in Qn have a capacity which is abnormally big, then the maximal flow (pn will be abnormally 
big too. The probability for these edges to have an abnormally large capacity is of order exp —Cn'^ 
for a constant C, because the number of edges in 0^ is C'n'^ for a constant C. On the opposite, if 
all the edges in a flat layer that separates T\ from in 0^ have abnormally small capacity, then 

will be abnormally small. Since the cardinality of such a set of edges is D'n'^~^ for a constant 
D', the probability of this event is of order exp— Dn"^"^ for a constant D. 

The second remark we would like to do is that the condition d{T^,T'^) > is relevant in Theorem 
[TOl First, without^ this condition, we cannot be sure that there exists a polyhedral set P as in 
the definition of <j)n, and thus the polyhedral approximation (see section [1| cannot be performed. 
Moreover, if ci(F^,F^) = 0, there exists a set of edges of constant cardinality (not depending on n) 
that contains paths from F^ to F^ through Qn for all n along the common boundary of F^ and F^, 
and so it may be sufficient for these edges to have a huge capacity to obtain that is abnormally 
big too. Thus, we cannot hope to obtain upper large deviations of volume order (see [26] for a 
counter-example). However, we do not know if this condition is essential for Theorem [1] to hold. 

3.2.1 Lower large deviations 

To prove Theorem [9l we have to study the probability 



for a positive e. The proof is divided in three steps. 

First step: We consider a set of edges £n that cuts F^ from F^ in of minimal capacity (so 
(pn = V{Sn)) and having the minimal number of edges among those cutsets. We see it as the (edge) 
boundary of a set En which is included in Q. Zhang's estimate of the number of edges in a minimal 
cutset (Theorem 1 in [28]) states that with high probability, the perimeter V{En,^) of En in $7 is 
smaller than a constant (5. Thus, En belongs to the set 



^ ^n< (0n - £)n' 



(1) 



Cf3 = {F cn\F cn, v{F, n)<p}. 



We endow with the topology L 



associated to the following distance d: 



diFi,F2) = £^(FiAF2), 
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where is the d-dimensional Lebesgue measure. For this topology, the set C/3 is compact. Thus, 
if we associate to each set F in a positive constant ep, and if we denote by V{F,eF) the neigh- 
bourhood of F of radius ep for the distance d defined above, the collection of these neighbourhoods 
is an open covering of Cp, and thus by compactness of Cjj we can extract a finite covering: 

N 



3Fi,...,Fn Cf3c\JV{Fi,eF,). 

i=l 

If we find an upper bound on the following probability: 

F [(An < {(Pn - £)n'^~^ and d(F„, F) < ep] (2) 
for each F in Cj3 and a corresponding ep, then we will obtain an upper bound on the probability 

Second step: We consider a fixed set F in C^, and we want to evaluate the probability 
So we suppose that is close to F for the distance d, we denote it by En ^ F to simplify the 
notations. We skip here all the problems of boundary conditions that arise in the proof of Theorem[9l 
we suppose that IniF) is equal to the integral of v along d*F n $7. 

We make a zoom along dF. Using the Vitali covering Theorem (Theorem [T2] in section [4]), we 
know that there exists a finite number of disjoint balls Bj = B{xj,rj) for j = 1, with Xj G dF 
such that dF is "almost flat" in each ball, and the part of dF that is missing in the covering has 
a very small area. We denote by vj the exterior normal unit vector of F at Xj (we suppose that it 
exists). Here "almost flat" means that 

(i) the capacity of dF inside Bj is very close to the capacity of the flat disc hyp(xj, vj) n Bj, i.e., 
very close to ad~irj~^i'{vj) ; 

(ii) F n Bj M BJ , where BJ is the lower half part of the ball Bj in the direction given by Vj-. 

BJ = {y G Bj\{y-Xj)-Vj < 0}. 

Thanks to property (i) and the fact that only a very small area of dF is missing in the covering, we 
know that 

Af 

In{F) is close to ^ aa-ir^-^H^j) ■ (3) 
i=i 

On the other hand, thanks to property (ii), we obtain that 

SnnBj ^ Fn Bj « BJ 

for the distance d. It means that in volume. En is very similar to BJ inside Bj, however there 
might exist some thin but long strands in Bj that belongs to En H {BjY or to En n BJ . We want 
to compare V{£n H Bj) with the maximal flow Tn{Dj,j) in a cylinder of basis Dj = disc{xj,rjvj) 
where r'j is a little bit smaller than rj, and 7 is a very small height, so that the cylinder is included 
in Bj and is almost flat. To make this comparison, we have to cut the above-mentioned strands 
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by adding edges to <f„. We do it very carefully, in order to control the number of edges we add, 
together with their capacity, and we obtain that 

ViSn DBj) < r„ {Dj , 7) + error , (4) 

where error is a corrective term that is very small. Combining ^ and Q, since Tq,{F) > (pn, we 
conclude that if < {cpn — e)n'^~^ and En ~ then there exists j G {1, AA} such that 

Tn{Dj,^) < {v{vj)-e/2)ad-irf-^n''-\ 



Third step: It remains to study the probability 

P[r„pj,7) < Hvj) - e/2)ad-irf'n^-^] . 

In fact it has already been done by Rossignol and Theret in [25]. It is easy to compare r„ (1)^,7) 
with a sum of maximal flows through cylinders whose bases are hyp err ect angles. Then, we can use 
directly Theorem 3.9 in [25] that states that the lower large deviations of these maximal flows below 
their limits are of surface order. 

3.2.2 Upper large deviations 

To prove Theorem [TOl we have to study the probability 



bn > {(t>n + £)n" 



(5) 



for a positive e. First of all, we can check that (j)Q is finite^In fact, we have to construct a polyhedral 
set P that satisfies all the conditions in the definition of (j)Q. This is done with the help of techniques 
very similar to some of those we will use in section |4] to complete our polyhedral approximation, so 
we will not explain these techniques here. The proof of theorem [10] is divided in three steps. 

First step: We consider a polyhedral set P as in the definition of 4>q such that IniP) is very 
close to this constant. We want to construct sets of edges near dP n Q that cut from F^ in 
iln- Because we took a discrete approximation of Q from the outside, we need to enlarge a little Q, 
because some flow might go from F^ to F^ using paths that lies partly in 0„ \ O. Thus we construct 
a set Q' which contains a small neighbourhood of Q (hence also Qn for all n large enough), which 
is transverse to dP, and which is small enough to ensure that Ifi'{P) is still very close to (pQ. To 
construct this set, we cover dQ with small cubes, by compactness we extract a finite subcover of 
dQ, and finally we add the cubes of the subcover to Q to obtain O'. We construct these cubes so 
that their boundaries are transverse to dP, and their diameters are uniformly smaller than a small 
constant, so that is included in a neighbourhood of Q as small as we need. Since dP is transverse 
to F, if we take this constant small enough, we can control Tf^-^idP n (O' \ 0)), and thus the 
difference between I^i{P) and Tn{P). 

Then we construct a family of Cn (where C > 0) disjoint sets of edges that cut F^ from F^ in 
and that lie near dP. We consider the neighbourhood P' of P inside at distance smaller 
than a tiny constant h, and we partition P' \ P into slabs M.'{k) of width of order 1/n, so we have 
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Cn such slabs which look like translates of dP n Vt' that are slightly deformed and thickened. We 
prove that each path from T\ to in 0^ must contain at least one edge that lies in the set Ai'{k) 
for each fc, i.e., each set M'{k) contains a cutset. Thus we have found a family of Cn disjoint cutsets. 

Second step: We almost cover dPr\ Q,' by a finite family of disjoint cylinders Bj, j £ J, whose 
bases are hyper rectangles of sidelength /, that are orthogonal to dP, of height bigger than h, and 
such that the part of dP which is missing in this covering is very small. Thus, we obtain that 

In'{P) is close to ^1/(1)^)/'^-^ (6) 

where Vj gives the direction towards which the cylinder Bj is tilted (it is the unit vector which is 
orthogonal to the face of dP that cuts Bj). 

We want to compare (j)n with the sum of the maximal flows (f){Bj,Vj). For each j, let Ej be a 
set of edges that cuts the top from the bottom of Bj. The set Uj^jEj does not cut F^ from F^ in 
r2„ in general, to create such a cutset we must add two sets of edges: 

(i) a set of edges that covers the part of dP n Q' that is missing in the covering by the cylinders 
Bj, 

(ii) a set of edges that glues together all the previous sets of edges (the sets Ej and the set 
described in (i)). 

In fact, we have already constructed Cn possible sets of edges as in (i): the edges that lie in 
M'{k) \ {Uj(zjBj) for k = l, Cn. We denote these sets by M{k). We can also find Cn (C > 0) 
disjoint sets of edges that can be the glue described in (ii), we denote these sets by W{1) for 
/ = 1, Cn. We do not provide a precise description of these sets. In fact, we can choose different 
sets because we provide the glue more or less in the interior of the cylinders Bj. Thus we obtain 
that 

VA; G V/ G {l,...,C"n} \J Ej U M{k) U W{1) cuts F^ from F^ in 0„ . 

J6J 



We obtain that 

"]cP{B„v,)+ min V{M{k)) + min V{W{1)) . (7) 

k=l,...,Cn 1=1,. ...Cn 



i6J 



Combining ^ and d?]), we see that if 0„ > {(pn +e)n , one of the following events must happen 



(a) 3j e J (l){Bj,Vj) > {u{vj) + e/2)l'^~^n'^-\ 

(b) VA;G {l,...,Cn} V{M{k)) > i]n'^-^ , 

(c) V/ G {l,...,Cn} V{W{1)) > r]n'^~\ 

where is a very small constant (depending on e and cpn). 

Third step: it consists in taking care of the probability that the events (a), (b) or (c) happen. 
The probability of (a) has already been studied in [26]: the upper large deviations of the variable 4> 
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in a cylinder above u are of volume order. The events (b) and (c) are of the same type, and their 
probability is of the form 



d-1 



m=l 



Dn 



(8) 



where (tm)meN is a family of i.i.d. variables of distribution function A, is a constant, ry is a 
very small constant and an'^~^ is the cardinality of the family of variables we consider. If a < 
r/E[ti]~^, and if the law A admits one exponential moment, the Cramer Theorem in M states that 
the probability ((8|) decays exponentially fast with n'^. Note the role of the optimization over Dn 
different probabilities to obtain the correct speed of decay. To complete the proof, it is enough to 
control the cardinality of the sets M{k) and W{1) for each k, I. This can been done, using the 
geometrical properties of dP (it is polyhedral and transverse to dil'). 



4 Polyhedral approximation : = 

We consider an open bounded domain Q in M'^. We denote its topological boundary by F = dQ. 
Let also F^, be two disjoint subsets of F. 

Hypothesis on U: We suppose that is a Lipschitz domain, i.e., its boundary F can be locally 
represented as the graph of a Lipschitz function defined on some open ball of M'^"^. Moreover there 
exists a finite number of oriented hypersurfaces Si, . . . , Sp of class which are transverse to each 
other and such that F is included in their union 5i U • • • U 5p. 

This hypothesis is automatically satisfied when is a bounded open set with a boundary or 
when $7 is a polyhedral domain. The Lipschitz condition can be expressed as follows: each point x of 
F = has a neighbourhood U such that UCiQ is represented by the inequality Xn < f{xi, • • • , 
in some cartesian coordinate system where / is a function satisfying a Lipschitz condition. Such 
domains are usually called Lipschitz domains in the literature. The boundary F of a Lipschitz 
domain is d — 1 rectifiable (in the terminology of Federer's book [H]), so that its Minkowski content 
is equal to H'^~^{T). In addition, a Lipschitz domain Q is admissible (in the terminology of Ziemer's 
book [29j) and in particular 7i'^"^{T \ 5*17) = 0. Moreover, each point of F is accessible from O 
through a rectifiable arc. 

Hypothesis on F^,F^: The sets F^, F^ are open subsets of F. The relative boundaries dr F^, 
dr F^ of F^, F^ in F have null 7i'^~^ measure. The distance between F^ and F^ is positive. 

We recall that the relative topology of F is the topology induced on F by the topology of M*^. Hence 
each of the sets F^,F^ is the intersection of F with an open set of W^. For F a subset of O having 
finite perimeter in O, the capacity of F is 

MF) = [ v{vF{y))dn''-\y)+ f v{vF{y))d'H''~\y)+ f y{v^^F{y)) dn<'-\y) . 
Jnnd*F Jr^nd'F Jr'^nd*{n-^F) 

o 

For all A C M'^, A is the closure of A, A its interior and A'^ = R"^ \ A. We will prove the following 
theorem: 
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Theorem 11. Let F be a subset of Q having finite perimeter. For any e > 0, there exists a 
polyhedral set P whose boundary dP is transverse to T and such that 

o ° 

T^CP, r2cM'^\P, c^iFA^Pnn)) < e, 



u{vp{x))dn'^-^{x) = IniP) < IniF)+e. 

d*Pnn 

First we notice that theorem [TT] implies that (pn = (j)n, and thus the convergence of cpn (see 
section [3^ . It is obvious since (pn < (see section [2^ . and theorem [TT] implies that (pn> (pn- 

The main difficulty of the proof of theorem [TT] is to handle properly the approximation close 
to r in order to push back inside all the interfaces. The essential tools of the proof are the 
Besicovitch differentiation theorem, the Vitali covering theorem and an approximation technique 
due to De Giorgi. Let us summarise the global strategy. 

Sketch of the proof: We fix 7 > 0. We cover d*Q up to a set of 7i'^~^ measure less than 7 by 
a finite collection of disjoint balls B{xi,ri), i G Ii U /2 U I3 U centered on T, whose radii are 
sufficiently small to ensure that the surface and volume estimates within the balls are controlled by 
the factor 7. The indices of /i correspond to balls centered on f^^*{Q.\ F), the indices of I2 to 
balls centered on F^ n d* F , the indices of Is to balls centered on (F \ F^) n d*F, the indices of 1^ 
to balls centered on (F \ F^) n 9* (17 \ F) (see figure [5]) . The remaining part of F is covered by a 

Balls indexed by I a 



Balls 
indexed 
by/i 



possible strands 
in dL 
for d > 3 




Balls indexed by 1^ 



Balls indexed by I a 



Balls 

indexed by 
h 



Figure 5: The balls indexed by Ij for i = 1, 5. 



finite collection of balls B{yj,Sj), j G Jq U Ji U J2. The indices of Ji correspond to balls covering 
the remaining part of Fi, the indices of J2 correspond to balls covering the remaining part of F2. 
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We choose e > sufficiently small, depending on 7 and on the previous families of balls and we 
approximate the set F by a smooth set L inside whose capacity and volume are at distance less 
than e from those of F. We build then two further family of balls: 

- B{xi,ri), i G Is, cover Vl n up to a set of l-f^"^ measure e. 

- B{yj, Sj), j £ J3, cover the remaining set Cl Pi dL \ Uje/5 B{xi, ri). 

Inside each ball B{xi,ri), i G Ii U /2 U I3 U /4 U Is, up to a small fraction, the interfaces are located 
on hypersurfaces and the radii of the balls are so small that these hypersurfaces are almost flat. 
Hence we can enclose the interfaces into small flat polyhedral cylinders D^, i G /i U /2 U Is U /4 U Is, 
and by aggregating adequately the cylinders to the set F or to its complement O \ F, we move 
these interfaces on the boundaries of these cylinders. The remaining interfaces are enclosed in the 
balls B{yj, sj), j G Jq U Ji U J2 U J3 and we approximate these balls from the outside by polyhedra. 

We have to define delicately the whole process, in order not to lose too much capacity, and to 
control the possible interaction between interfaces close to T and interfaces in Q. The presence of 
boundary conditions creates a substantial additional difficulty compared to the polyhedral approx- 
imation performed in [9]. Indeed, the most difficult interfaces to handle are those corresponding to 
Di, i e I3U I4. We first choose the balls B{xi,ri), i e Ii U I2 U U I4, corresponding to 7. We 
cover the remaining portion of F with the balls B{yj,Sj), j G Jq U Ji U J2. At this point we can 
already in principle define the cylinders Dj, i G /i U/2- Then we choose e small enough, depending 
on 7 and the balls B{xi, rj), i G Ii U I2 U /a U I^, to ensure that the perturbation of volume e caused 
when smoothing the set F inside O will not alter significantly the situation inside the balls B{xi,ri), 
i G Is U /4 . Then we move inside Q and we build the cylinders Di, i £ I^. Then we come back to 
the boundary and we build the cylinders Dj, i G Is U /4. We cover the remaining interfaces in O 
by the balls B{yj,Sj), j G J3. Finally we aggregate successively each flat polyhedral cylinder Di to 
the set L or to its complement. 

Preparation of the proof. Let us consider a subset F of^l having finite perimeter. Let 7 belong 
to ]0, 1/16[. We start by handling the boundary F, for which we make locally flat approximations 
controlled by the factor 7. By hypothesis, there exists a finite number of oriented hypersurfaces 
5*1, . . . , 5p of class such that F is included in their union 5i U • • • U Sp. In particular, we have 

T\d*n C S = (j SkDSi. 

l<k<l<p 

Since the hypersurfaces 5i, . . . , 5^ are transverse to each other, this implies that H'^~^{S) = 0. 

• Continuity of the normal vectors. The hypersurfaces Si,...,Sp being and the set F 
compact, the maps 2; G F vs^{x), 1 < k < p (where vsf.{x) is the unit normal vector to Sk at x) 
are uniformly continuous: 

V(5> 37?>0 VA: G {!,... ,p} Vx,yG5fcnF \x - y\2 < t] ^ \vs,{x) - vs,,{y)\^ < 5 . 

Let T]* be associated to 5 = 1 by this property. We will use also a more refined property. 

• Localisation of the interfaces. We first prove a geometric lemma: 

Lemma 1. Let F be an hypersurface (that is a suhmanifold ofW^ of codimension 1) and let K 
he a compact subset ofT. There exists a positive M = M{T,K) such that: 

Ve > 3r>0 \/x,y e K \x - y\2 < r =^ d2{y , tan{T , x)) < M e \x - y\2 ■ 

(tan{T,x) is the tangent hyperplane ofT at x). 
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Proof. By a standard compactness argument, it is enough to prove the following local property: 

Vxer 3M{x)>0 Ve>0 3r(x,e)>0 ^y, z e F n B{x,r{x,e)) 

d2{y, tan(r, z)) < M{x) e\y — z\2 ■ 

o 

Indeed, if this property holds, we cover K by the open balls B{x,r{x,e)/2), x G K, we extract a 
finite subcover B{xi,r{xi, e)/2), 1 <i < k, and we set 

M = max{ M{xi) : 1 < i < A: } , r = min{ r{xi, e)/2 : 1 < i < k} . 

Let now y,z belong to K with \y — z\2 < r. Let i be such that y belongs to B{xi,r{xi,e)/2). Since 
r < r{xi,e)/2, then both y,z belong to the ball B{xi,r{xi,£)) and it follows that 

c?2(y,tan(r,z)) < M{xi)£\y-z\2 < Me\y-z\2. 

We turn now to the proof of the above local property. Since T is an hypersurface, for any x 
in r there exists a neighbourhood F of a; in W^, a diffeomorphism / : y i— > M'^ of class and a 
{d — 1) dimensional vector space Z of such that Z n f{V) = f{T n V) (see for instance |14| . 
3.1.19). Let ^ be a compact neighbourhood of x included in V. Since / is a diffeomorphism, the 
maps y £ A df{y) £ End(M'^), u G f{A) df~^{u) G End(M'^) are continuous. Therefore they 
are bounded: 

3M>0 yyGA \\dfiy)\\<M, yu€f{A) \\df-\u)\\ < M 

(here ||(i/(x)|| = sup{ \df{x){y)\2 : \y\2 < 1 } is the standard operator norm in End(M'^)). Since f{A) 
is compact, the differential map df~^ is uniformly continuous on f{A): 

Ve> 36 >0 yu,vef{A) \u - v\2 < S ||d/~^(n) - 4f"^(t;)|| < e . 

Let e be positive and let 6 be associated to e as above. Let p be positive and small enough so 
that p < S/2 and B{f{x),p) C f{A) (since / is a diffeomorphism, f{A) is a neighbourhood of 
f{x)). Let r be such that < r < p/M and B{x,r) C A. We claim that M associated to x and r 
associated to e,x answer the problem. Let y,z belong to F n B{x,r). Since [y,z] C B{x,r) C A, 
and \ \df{C)\\ <MonA, then 

\f{y) - f{x)\2 < M\y - x\2 < Mr <p, \f{z) - f{x)\2 <p, 
\f{y) - f{z)\2 < 6 , \f{y) - f{z)\2 < M\y - z\2 . 

We apply next a classical lemma of differential calculus (see [20], I, 4, Corollary 2) to the map 
and the interval [f{z),f{y)] (which is included in B{f{x),p) C f{A)) and the point f{z): 

\y-Z-dr\f{z)){f{y)-f{z))\2 < 

\f{y) - /(z)|2sup{ ||#-i(C) - dr\f{z))\\ : C G [f{z)J{y)] } . 

The right-hand member is less than M\y — z\2e. Since z + df'^{f{z)){f{y) — f{z)) belongs to 
tan(r,2;), we are done. □ 

We come back to our case. Let k G {l,...,p}. The set 5^ n F is a compact subset of the 
hypersurface Sk- Applying lemma [H we get: 

3Mk V^o > 3 ?7fc > V.X, y G Sfe n T |x - yjs < 7?fc ^ d2{y, tan(S'fc, x)) < MkSo\x - y\2 • 
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Let Mq = maxi<k<p Mk and let 60 in ]0, l/2[ be such that M060 < 7- For each /c in { 1, . . . ,p }, let 
r/fc be associated to 60 as in the above property and let 

r?o = min ( min 77*, ■^dist(r\r2; 

• Covering of F by transverse cubes. We build a family of cubes Q{x,r), indexed by x S F 
and r £]0,rr[ such that Q{x,r) is a cube centered at x of side length r which is transverse to F. 
For X eR'^ and k e {I, . . . ,p}, let Pk{x) be a point of 5"^ n F such that 

\x - Pk{x)\2 = inf {\x - y\2 : y £ Sknr} . 

Such a point exists since 5^ n F is compact. We define then for k £ {1, . . . ,p} 

Vx G M"' Vkix) = vs^{pk{x)). 

We define also 

dr = inf max min (\e — Vi\2,\ — e — Vi\2) 
vi,...,vp&S''-'^ beBd I < k < r 

e£b 

where Bd is the collection of the orthonormal basis of and S'^~^ is the unit sphere of M*^. Let rj 
be associated to dr/4: as in the above continuity property. We set 

T] 

rr = — ; • 
2d 

Let X E F. By the definition of dr, there exists an orthonormal basis bx of M.'^ such that 

\/e£bx VA; G { 1, . . . ,p} min {\e - Vk{x)\2,\ - e - Vk{x)\2) > ^ ■ 

Let Q{x,r) be the cube centered at x of sidelength r whose sides are parallel to the vectors of bx- 
We claim that Q{x,r) is transverse to F for r < rp. Indeed, let y G Q{x,r) n F. Suppose that 
y G Sk for some A; G { 1, . . . so that Vk{y) = vSkiv) \x — Pk{x)\2 < dr^. In particular, we 
have \y - Pk{x)\2 < 2drr < t] and \vs^{y) - ffc(x)|2 < dr/i. For e G bx, 

Y < \e-Vk{x)\2 < \e- vs^{y)\2 + \vs,,{y) - Vk{x)\2 

whence 

I /■ \\ ^1" 

\e-vs^{y)\2 ^ y - ^ = X' 



This is also true for — e, therefore the faces of the cube Q{x,r) are transverse to Sk- 
• Vitali covering Theorem for 7i'^~^. A collection of sets U is called a Vitali class for a Borel set 
E of M'^ if for each x £ E and 6 > 0, there exists a set U (zlA containing x such that < diam U < 6, 
where diamf/ is the diameter of the set U. We now recall the Vitali covering Theorem for 7i'^~^ 
(see for instance [13], Theorem 1.10), since it will be useful during the proof: 
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Theorem 12. Let E be a Ti'^ ^ measurable subset of and Li be a Vitali class of closed sets for 

E. Then we may select a (countable) disjoint sequence {Ui)i^i from Li such that 

either ^(diam C/j)"'-^ = +00 or \ Uje/C/j) = 0. 

is/ 

jd-ii 



IfTi (E) < 00, then given e > 0, we may also require that 

2d 



Start of the main argument. We first handle the interfaces along T. Let 7^(r) be the set of the 

points 2; of r \ such that 



lim (ad/)-^£'^(S(x,r) \!^) = 1/2, 

r— >0 



lim {ad-ir''~T'n''^\B{x,r)nr) = 1. 
Let TZ{n \ F) be the set of the points x belonging to d*{n \ F) n 7^(^) such that 



lim {ad^ir'^-Y^n'^~\B{x,r)nd*{n\F)) = 1 
1 — >o 

d\-l rdi 



lim {adr'')-'C''{B{x, r) n (0 \ F)) = 1/2 , 
lim^ [ad^ir^-^r^ [ iy{vn~^F{y)) dU^-^y) = u^x)) . 



r-->0 



B(x,r)r\d*(n-^F) 



Let 'Tl{F) be the set of the points x belonging to (9*F n 7^(r) such that 



lim {ad-ir'^~^y^'H'^^^{B{x,r)r\d''F) = 1 
1 — >o 



lim {adr'')-'C''{B{x,r)r\F) = 1/2, 



lim iad-ird-^' / ,,{vF{y)) dn^-'iy) = i^ivn{x)) . 

^'^^ JB{x,r)r\d*F 

Thanks to the hypothesis on T and the structure of the sets of finite perimeter (see either Lemma 1, 
section 5.8 of [12], Lemma 5.9.5 in [29] or Theorem 3.61 of [I]), we have 

n'^-^ (r \ (7^(F) u n{n \ f))) = o . 

For X in 7^(r), there exists a positive ro(x,7) such that, for any r < ro(rE,7), 

\C'^iB{x,r)\n)-adr'^/2\ < adr\ 
\n'^-\B{x,r)nT)-ad^ir'^-^\ < jad-ir'^-^. 

For X in 7^(17 \ F), there exists a positive r{x,^) < ro(x,7) such that, for any r < r(x,7), 
\n''-\B{x,r)nd*{n-^F))-ad-ir''-'\ < ^aa-^-^ 

[ad-ir^-^r^ I y{vn^F{y))d-H''-\y)-v{vn{x)) < 

J B{x,r)r\d*{n-^F) 
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For X in TZ{F), there exists a positive r{x,j) < ro(x,7) such that, for any r < r{x,j), 

\n'^-\Bix,r) nd*F) - ad^ir'^~'\ < j ad-ir'^^^ 
\£'^{B{x,r)nF)-ay/2\ < ^ , 



JB{x,r)r\d*F 



< 7 



Let us define the sets 



r^* = n n{n \ f) , r^* = r^n n{F) , 
p3* = (r \ fa) n TZ{F) , r^* = (r \ fi) n n{n \ f) . 



The family of balls 



B{x,r) , xGr^*ur^*, r < min (r(x, 7), 7, ?7o, ^dist(2;, 5)) , 

1 1 — 

B{x,r), xeV^*, r < min (r(a;,7),7,r?o, -dist(rE, 5), -dist(2;,r2)) , 

B{x,r), xeT^*, r < min (r(j;,7),7,?7o, ^dist(rE,5), ^dist(j;,ri)) 

is a Vitah relation for T^* U T^* u T^* U T^* . Recall that S is the set of the points belonging to 
two or more of the hypersurfaces Si,...,Sp and since S is disjoint from pi*^ ^2* -p3* p4*^ then 
dist(a;, 5) > for x E T^* U T^* U T^* U T'^* . By the standard Vitali covering Theorem (see theorem 
[T2l) . we may select a finite or countable collection of disjoint balls B{xi,ri), i £ I, such that: for 
ieI,Xie ri* ur^* ur^* ur^*, < min(r(xi,7),7,r/o, idist(xj,5)) and 

either H'^^^ (t \\J B{xi,ri)^ = or ^rf-^ = oo. 

iGl iel 

Because for each i in /, is smaller than r(xj,7), 

arf_i(l-7)5]rf-i < H'^-Ht) < 00 

and therefore the first case occurs, so that we may select four finite subsets Ii, I2, 13, h of / such 
that 

VA: G {!,... ,4} Vi G 4 oci e T''* , 
W^-i(rx (J IJ B{xi,u)) < 7. 

l<A:<4ie/fe 

Let i belong to /i U /a U /a U 14. We have 

n''^\r n n) \ r,(i - 2^))) = ^^-^(r n B{xi, r^)) - n''-\r n ri(i - 2^))) 

< (1 + 7)«d-irf"i _ (1 _ _ 2^)'^-! 

= arf_irf"i(l + 7 - (1 - 7)(1 - 2^7)'"') 



23 



4 POLYHEDRAL APPROXIMATION 



Hence 

^ n''-\rnB{xi,ri) \ Bix^,n{i - 2^))) 

ie/iU/2U/3U/4 

ie/iU/2U/3U/4 

and 

We have a finite number of disjoint closed balls rj(l — 2^/7)), i G /i U/2 U/3U/4. By increasing 
slightly all the radii rj, we can keep the balls disjoint, ensure that each radius satisfies the same 
strict inequalities for i in /i U /2 U /a U /4, and get the inequality 

The above set is a compact subset of T. For k = 1,2, we define 

Rk = Tk\ U B{x„u{l-2^)). 

The sets Ri and R2 are compact and their 7i'^~^ measure is less than 27 + 4d^7i'^~^{T) (recall 
that drT^ and drl"^ have a null 7i'^~^ measure). For A; = 1,2, by the definition of the Hausdorff 
measure 7i'^~^, there exists a collection of balls B{yj, sj), j G such that: 

yj £ Jk < < min (r/o, y) , B{yj, sj) n Rk ^ , 
^fc C IJ B{yj,Sj). 

By compactness of Ri and R2, the sets Ji and J2 can be chosen to be finite. It remains to cover 

i?o = r \ U B{xi,uil - 2^7)) \ U • 
iG/iU/2U/3U/4 ieJiUJ2 

The set Rq is a closed subset of F which is at a positive distance from F^ and F^. There exists a 
collection of balls B{yj,Sj), j £ Jq such that: 

VjGJo 0<Sj <min(r/o,Y,^dist(i?o,r^UF2)), Biyj, sj) D Rq ^ , 

ad-is'f' < 37 + M^n'^'H^) , 

J6J0 

Ro C \J B{yj,Sj). 
jeJo 
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Now the collection of balls 

B{x„ri{l - 2^)), iGhUhUhU h, B{yj, sj), j G Jq U Ji U J2 

covers completely T. We will next replace these balls by polyhedra. For j € Jq U Ji U J2, let Xj 
belong to B{yj,Sj) H T and let Qj be the cube Q{xj,4:Sj). For i in Ii U I2 D Is L) 1^, the point Xi 
belongs to exactly one hypersurface among Si, . . . ,Sp, which we denote by Ss{i)- In particular F 
admits a normal vector vn{xi) at Xi in the classical sense. For each i in /i U /2 U /s U I4, let Pi be a 
convex open polygon inside the hyperplane hyp{xi,vn{xi)) such that 

disc{xi,ri{l - 2^),vn{xi)) C Pi C disc(xi,ri(l - ^/^),vn{xi)) , 

\n''-\n) - ad-irf-i(l - V7)'~'l < Soad-irt\l - ^f-^ . 
Thanks to the choices of the radius rj and the constants Mq, r/o, we have then 

rnB(xi,r,(l-2^)) C %)n5(xi,ri(l-2^)) C cyl(Pi,27ri), 
Tr^B{xi,ri) C Ss^i)r^ B{xi,r-^ C cy\{(Msc{xi,ri,vn{xi)), Midori) , 
Vx £ B{xi,ri) nT \vn{x) - vn{xi)\2 < 1. 

The choice of 5o guarantees that Mo6o{l + So)ri < 2^ri. Let t be such that 

Modoil + do)ri < t < ^ri . 

We have 

-tvn{xi) + Pi cnn B{xi, r,) , F n {-tvn{x^) + Pi) = . 

In particular, the set F can intersect the cylinder cyl(Pj,t) only along its lateral sides, which are 
parallel to vn{xi). Let x belong to F n 9cyl(Pj,t). Then 

\vcy\(Pa)i^) - > \vcyi{P,,t){x) - vn{xi)\2 - \vn{xi) - vn{x)\2 > V2-I. 

o 

Therefore the cylinder cy\{Pi,t) is transverse to F. We will replace the ball B{xi,ri{l — 2.^/7)) by the 
cylinder cyl(Pi, tj), for a carefully chosen value of tj in the interval [Mo5o(l + So)ri, ■y/7rj[. However, 
we must delay the choices of the values ti, i £ I3U I4 until we have modified the set F inside Q. We 
deal next with the interfaces inside and we make an approximation of F controlled by a factor 
e. We choose e sufficiently small compared to 7 so that, when we perturb the set F by a volume e, 
the resulting effect close to F is still of order 7. Let e be such that < e < 7 and 

e < 'jad min rf . 

We use next a classical approximation result: there exists a relatively closed subset L of O having 
finite perimeter such that $7 n dL is an hypersurface of class C°° and 



C'iFAL) <e, [ v{vF{y))dn''~\y)- [ v{vL{y)) dn''-\y) 

Jnnd*F JnndL 



< e. 
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In the case where u is constant, this result is stated in Lemma 4.4 of [23]. In the non constant case, 
the argument should be slightly modified, as explained in the proof of proposition 14.8 of [9], where 
the approximation is performed in M'^ instead of Q. When working inside Q, the extra difficulty is 
to deal with regions close to the boundary (see the proof of Proposition 4.3 of [23]). For r > 0, we 
define 

dLr = {x £ dL : d{x, T) > r } . 
By continuity of the measure 'H'^~^\dL, there exists r* > such that 

n'^~^{nndL\dL2r*) < e. 

We apply lemma [U to the set dLj.* and the hypersurface VL n dL: 

3M > V(^> 3 7?>0 yx,y£dLr-' \x - y\2 < r] ^ d2{y,tan{dL, x)) < M5\x - y\2 . 

For a point x belonging to dLr*, the tangent hyperplane of Pi dL at x is precisely hyp(x, 

Let M be as above. We can assume that M > 1. Let 6 in ]0, (^o[ be such that 2SM < e. Let r\ be 

associated to b as in the above property. For x G dL2r* , 

lim {ad-ir'^-^Y^n'^~^[B{x,r)^dL) = 1, 
hm^ (a,_ir'^-i)-i / vivL^y)) dn''-\y) = v{vl{x)) . 



B{x,r)ndL 



For any x in dL2r* , there exists a positive r(a;,e) such that, for any r < r{x,e), 

\n'^~^{B{x,r)ndL) -ad-ir'^~^\ < ead_i/"\ 
'ad-ir^-Y' [ u{vL{y))dn''-\y) - u{vl{x)) < e . 

J B{x,r)r\dL 

The family of balls B{x,r), x G dL2r*, r < min(r*, r/o, e), e, 77), is a Vitali relation for dL2r*- 
By the standard Vitali covering Theorem, we may select a finite or countable collection of disjoint 
balls B{xi,ri), i £ /', such that: for any i in /', Xi S dL2r* , 

ri < mm{r*,r]o,r{xi,e),e,r]) 

and 

either rL'^^^(dL2r* \ |J B{xuri)^ =0 or ^r'^^'^ = 00. 

iei' iei' 

Because for each i in /', is smaller than r{xi,e), 

ad-i{l-e)^rf~^ < W^-^nndL) < 00 
ie/' 

and therefore the first case occurs, so that we may select a finite subset of /' such that 

n''^^{dL2r* \ IJ B{xi,n)) < e. 
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We have a finite number of disjoint closed balls B{xi,ri), i G I5. By increasing slightly all the radii 
Tj, we can keep the balls disjoint, each strictly smaller than min(r*, r/o, ^(xj, e), e, ry) for i in Is, 
and get the stronger inequality 

For each i in /s, let Pi be a convex open polygon inside the hyperplane hyp(xj, ^^(xi)) such that 

disc{xi, ri,VL{xi)) C Pi C disc(xj,ri(l + 5),VL{xi)) , 

\n''-\P.)-ad-irt^\<6ad-irt^. 

We set ijj = M5{1 + 5) (hence ■0 < £ < !)• Let i belong to Is. Let Di be the cylinder 

Di = cy\{Pi,M5{l+5)ri) 

of basis Pi and height 2'i/'rj. The point Xj belongs to dL2r*, the radius is smaller than and r* , 
so that 

\/x e dLn B{xi,ri) d2{x,h.yp{xi,VL{xi))) < M6\x - Xi\2 , 

whence 

o 

dLnB{xi,ri) C cyl (disc(xi,rj,-!;L(2;i)),M(5ri) C A- 

We will approximate F by L inside O and we will push the interfaces rid*{Q\ F) and F^ n d*F 
into ri. We next handle the regions close to F inside the family of balls B{xi,ri), i G /1U/2U/3U/4. 
We will modify adequately the set F to ensure that no significant interface is created within these 
balls. Our technique consists in building a small flat cylinder centered on F which we add (for 
indices in Ii U I3) or remove (for indices in I2 U I4) to the set F. We have to design carefully this 
operation in order not to create any significant additional interface. This is the place where we tie 
together the covering of the boundary and the inner approximation. Recall that we already chose 
a family of polygons Pj, i G IiU I2D I3U I4. For i £ IiU I2, we simply define Di to be the cylinder 

Di = cyl{P^,Mo6o{l + So)ri), 

see figure [6l The construction of the cylinders associated to the indices i £ I3 U I4 is more com- 
plicated. Our technique consists in choosing carefully the height U of the cylinders cyl(i-'j,tj) for 
i £ I3U I4. We examine separately the indices in I3 and I4,. 

• Balls indexed by I3. Let i belong to I^. Because of the condition imposed on e, we have 
\C'^{B{x„ri)nL)-adrf/2\ < j a^rf + e < 2-f a^rf . 

Since in addition 

\£''{B{xi,r,) \ f]) - adrf/2\ < 7 a^rf , 

it follows that 
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Thanks to the choice of the polygon Pi, we have then 

/ n'^-\{-tvn{xi) + P^) \L)dt < C\B{x^,n) n \ L)) < ^^advf . 

The condition on 7 yields in particular y/^ — 27 > ^/2. Hence there exists ti s]27rj, ^ri[ such 
that 

n''-\{-Uvn{xi) + P^) \ L) < e^a^rf^i . 
Let Di be the cylinder Di = cyl(Pj,tj). 

• Balls indexed by I^. Let i belong to I4. Because of the condition imposed on e, we have 
\C\B{xi,n)r^{^\L))-adrf/2\ < ja^rf + e < 2^ a^rf . 

Since in addition 

\£''{B{xi,r,)^n)-adrf/2\ < j a^rf , 

it follows that 

C\Bixi,r,)nL) < S^aarf. 
Thanks to the choice of the polygon Pi, we have then 

f n'^-^{{-tvn{xi)+Pi)r\L)dt < C'^{B{xi,n)r\L) < Sjaarf. 

The condition on 7 yields in particular ^ — 27 > ■y/7/2. Hence there exists ti €]27rj, ^ri[ such 
that 

n''~\{-tivn{x^) + Pi) n L) < 6^adrf-i . 

Let Di be the cylinder Di = cyl(Pj,ii) (see figure [7|. We have now built the whole family of 
cylinders Dj, i G Ii U /2 U Is U /4 U Is. Moreover, the sets 

A, i G/1U/2U/3U/4, B{yj,Sj), j G J0UJ1UJ2, 
cover completely T. It remains now to cover the region 

i?3 = ^ndL \ U ^ U ^ivj^^j)- 

iehuhulsUliUls jeJoUJiUJ2 
Since R3 does not intersect F, the distance 

P = ^dist(r,ii3) 

is positive and also R3 is compact. Prom the preceding inequalities, we deduce that 

T-c^^^Rs) < n'^-^{nr\dL\dL2r*) + n'^^^(dL2r* \ U A 

< e + n'^~^(dL2r' \ U ^(xi,ri)) < 2e. 
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Figure 7: The cylinder Di for i € 14^. 
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By the definition of the Hausdorff measure H ^, there exists a collection of balls B{yj, Sj), j € J3, 
such that: 

Vj G J3 0<Sj < p, B{yj,Sj) n i?3 / , 

By compactness, we might assume in addition that J3 is finite. For j € J3, let Xj belong to 
B{yj, Sj) n i?3 and let Qj be the cube Q{xj,Asj). We set 

P = ({nnL)u [j A u U Qj^ \ U A \ U Qj- 

^ i&huhuh jeJi ^ iehuh jeJoUJ2UJ3 

The sets 4', J ^ -^0 U Ji U J2 U J3, A, i G /i U /2 U 13 U /4 U h cover dL U T, therefore 

dP c U aAu U ^^J ' 

iG/lU/2U/3U/4U/5 igJoUJlUJ2UJ3 

thus P is polyhedral and dP is transverse to T. Since the sets 

A , i G /i U /3 , 4 ' J e -^1 
cover completely F^, while the sets 

Di, iehuhuh, Qj, j G Jo U J2 U J3 
do not intersect F^, then F^ is included in the interior of P. Similarly, the sets 

A , iehuh, Qj, j e J2 

2 

cover completely F , while the sets 

Di, iehuhuh, Qj, j G Jo U Ji U J3 
2 2 

do not intersect F , thus F is included in the interior of the complement of P. We next check that 
the set P n ri approximates the initial set F with respect to the volume. We have 

{Pnn)AF C {LAF)U (J Au U Qj 

iehuhuhuhuls ieJoUJiUJ2UJ3 

whence 

C'^i{Pnn)AF) < e+ 

J2 2Qrf_irf-i(l + 5o)^/7^* + 5^2arf_irf-^(l + 5)^ri+ J2 M'^Sjf . 
iehuhuhuh i&h jG J0UJ1UJ2UJ3 
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Yet each rj is smaller than 7, 

Y^aa^irf-' < 2n''-\Qr\dL) < {u^^^n''~\d* F n U) + e) , 
ad-is'^-^ < 3(37 + Ad^n''~\r)) + 3e , 

jeJoUJiUJ2UJ3 

so that 

^min 

+ 3 . 2^^(37 + Ad^'H''~\Y) + e) . 

We estimate next the capacity of P. To do this, we examine the intersection of dP n with each 
polyhedral cylinder. For i G Ji U /2, we use the obvious inclusion 

Pr\^r\dDi c nndDi. 

For i G /3U/4U/5, the sets 9Pnr2n9Dj require more attention. We consider separately the indices 
of Is, I 4 and I5. 

• Cylinders indexed by I3. Let i in I3. We have 

nndPndDi c nn (aA \ u (J sg^- . 

jreJoUJiUJ2UJ3 

Yet, thanks to the construction of the cylinder A, 

n'^-^nndDi^L) < n'^-\{-tivn{xi) + Pi) ^i) +n'^-\dPi)2^ri 

• Cylinders indexed by 1^. Let i in 14. We have 

nndPndDi c sin (aA n l) u |J 9Qj . 

ieJoUJiUJ2UJ3 

Yet, thanks to the construction of the cylinder A, 

n'^~\nndDinL) < n'^~\{-uvnixi) + Pi)nL) + n'^~\dPi)2^ri 

< d^adrt^ + 2ad^2rt^2^ri < 6^(ad + ad^2)rt^ ■ 

• Cylinders indexed by I5. Let i in Is. We set 

Gj = disc (xj - ipnvLixi), ip'^ri,VL{xi)) . 

We claim that the set Gi is included in the interior of L. Indeed, Gi C B{xi,ri) n 9A, yet 
9-L n B{xi,ri) C A, therefore Gi does not intersect 5L. Since VL{xi) is the exterior normal vector 
to L at Xi, then Gi is included in L. The definition of the set P implies that 

dPnGi C \J dQj, 
ieJoUJiUJ2UJ3 
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whence 



QndPndDi c {dDi\Gi)U |J dQj . 

jeJoUJiUJ2UJ3 



Yet 



W-^-i (aA \ {Pi + i^TiVLixi)) \ Gi) < 2ad-2rt^2iPri + a^-irf-^ {l + 5 - {1 - ^2)('i-i)/2) 

< a,.^rt'U^^P + I + 6 - {1 - . 

V aa-i ) 



Finally, we conclude that 



n 9P c U (5^ n dDi) u U (rj n A \ u J (0 n aA n L) 
U IJ (9A \ G,) u U . 

ie/s j6JoUJiUJ2UJ3 



Therefore 



~l~ ^max 

ie/4 

ieJoUJiUJ2UJ3 
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We use now the various estimates obtained in the course of the approximation. We get 

iehuh 

iei'3U/4 



+ Y[ad-irt\l + S)u{vL{xi)) 

+ -n.ax«.-irf-^ + 1 + 5 - (1 - 

+ 1^T. I v{vL{y))dH''-\y) 



Ea_i / ad-2 ^ . „ ^ad + 0!d^2 , . ad-2 . 
t'lnaxttd-l'^i —J 57 + h 4 V 

+ 1 + 5 - (1 - ^2)(rf-l)/2^ + j.^^^2'^-i3(37 + id^n''-\T) + e) 



< / Hvn{y))dn''~\y)+ [ u{vn{y))dn^-\y) 

+ / i^ivLiy))dn''~Hy) 

JnndL 

+ 2{n'-HT) + n'~\n n dL))u^J^Bj + 6^7^^^^ + 

+ 1 + 5 _ (1 _ ^2)(d-i)/2^ ^ ^.„,ax(2''-^3(37 + 4(i^W'^-i(r)) + 3e 

1 + 60 

7 



1-7 



+ ^(r)H — U'max 57 + 6^ hfe + 4- 



-ad-1 ad-1 Od-i 

+ i^max (2^-^3(37 + Ad^H'^-Hr)) + 3e^ 

where we have used the inequahty ip < e in the last step. We have also use the inclusions 

ViG/i B{xi,ri)nd*{n\F) c r^na*(0\F), 

Since 5oi^)7i£ can be chosen arbitrarily small, we have obtained the desired approximation. □ 
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5 Positivity of (jjn 

We suppose that 

/ xdA{x) < oo , (9) 

^[0,+oo[ 

We will prove that (pfi > if and only if A(0) < 1 — Pc{d). In fact we know that if the condition ^ 
is satisfied, 

A{0)< I -pc{d) <^ 3v,u{v)>0 \/v,u{v)>0. 

Thus, the implication ^ 

A{0)>l-p,{d) =^ cbn = 

is trivial. We suppose that A(0) < 1 — pc{d). Since u satisfies the weak triangle inequality, the 
function v i— > ^'{v) is continuous, and so as soon as A(0) < 1 — Pc{d) and ([9]) is satisfied, we have 

f mill = min u > . 

If P is a polyhedral set, then n'^-^{{dP nn)\ {d*P n ft)) = 0. We then obtain that 
(h > J^min X mf{n'^^^{S nn)\S hypersurface that cuts from in H, d{S, U T^) > 0} . 

We recall that the hypersurface S cuts from in O if 5 intersects any continuous path from a 
point in to a point in F^ that is included in 0. We consider such a hypersurface S C M^, and we 
want to bound from below the quantity H'^~^{S H Q) independently on S. 

The idea of the proof is the following. We consider a path from F^ to in Q. We construct a 
tubular neighbourhood of this path of diameter depending only on the domain and not on the path 
itself that lies in except at its endpoints. Then we prove that it is not very deformed compared 
to a straight tube. Since S has to cut this tube, we obtain the desired lower bound 7i'^~^{S H Q). 

For i = 1, 2, we can find Xi in F* and > such that F n B{xi, ri) C F* and F n -B(xj, ri) is 
a hypersurface. We denote by vn{xi) the exterior normal unit vector to i7 at Xj, and by Tii{xi) 
the hyperplane tangent to F at Xj. Since F is of class in a neighbourhood of Xi and O is a 
Lipschitz domain, applying lemma [H we know that for all ^ > 0, there exists e > depending on 
(0, F, F^, F^, xi, X2) such that for i = 1, 2 we have 

VLf\B{xi, 2e) is connected , 
Tr\B{x„2e) C V2{Tn{xi),2esine)r\B{xi,2e), 
Vr\B{xi,2e) C T\ 

We fix 6 small enough to have 2esin0 < e/2. We define 

Ai = TQ,{xi)r\B{xi,e) and A = cyl(^i,e), 

and then 

where Di is the interior of Di for z = 1,2. We define 

Xi = {z e bi \ XiZ ■ vn{xi) > e/2} c O . 



35 



5 POSITIVITY OF 4>n 



Then Xi C Q \ ^l. Each path r from a point yi € Xi to a point y2 £ X2 contains a path r' from a 
point y'l G to a point y'2 G such that r' C 0,, thus S intersects r. We consider the set 

Vi = {z£Xi\d2iz,dX,)>e/8}. 

Let yi G Vi, ^2 £ V2 such that d2{yi, dXi) > e/A for z = 1, 2. Since is obviously connected by arc, 
there exists a path r from yi to 1)2 in VL. The path f is compact and O is open, so 5 = d2{r, dO,) > 0. 
We thus can find a path r included in V2(r, min((5/2, e/8)) which is a C°° submanifold of of 
dimension 1 and which has one endpoint, denoted by yi, in Vi, and the other one, denoted by y2, 
in V2. ^ ^ 

As we explained previously, d2{r,d^l) > 0, so there exists a positive rji such that V2(r, r/i) C ^l. 
We can suppose that rji < e/16, to obtain that B{yi, rji) C Xi for i = 1,2. For all z in r we denote 
by Nr{z) the hyperplane orthogonal to r at z, and by Nr{z) the subset of Nr{z) composed of the 
points of Nr{z) that are at distance smaller than or equal to rj of z. The tubular neighbourhood 
of r of radius rj, denoted by tub(r, r/), is the set of all the points z in such that there exists a 
geodesic of length smaller than or equal to 77 from z that meets r orthogonally, i.e., 

tub(r,r?) = IJ^r(^)' 

(see for example [l7]). We have a picture of this tubular neighbourhood on figure [H Since r is a 
compact C°° submanifold of M.'^ which is complete, there exists a 7^2 > small enough such that for 
all T] < rj2, the tubular neighbourhood of r of diameter r/ is well defined by a C°°-diffeomorphism 
(see for example [3], Theorem 2.7.12, or [l7]), i.e., there exists a C°°-diffeomorphism V' from 

Nr'^ = {{z, v), z £r,v £ N^iz)} 

to tub(r, r/). We choose a positive t] smaller than min(r/i, 772)- We stress the fact that this rj depends 
on (0,F,F\f2) but not on S. 

Let {I,h) be a parametrisation of class C°° of r, i.e., I = [a,b] is a closed interval of M, 
/i : / — > r is a C°°-diffeomorphism which is an immersion. Let z be in r, and = h~^{z) G I. 
The vector h'{uz) is tangent to r at 2, and there exists some vectors {e2{z), ...,ed{z)) such that 
{h'{uz),e2{z), ...,ed{z)) is a direct basis of W^. There exists a neighbourhood Uz of Uz in / such 
that for all u G Uz, {h'{u),e2{z), ...,ed{z)) is still a basis of M"^, since h' is continuous. Indeed 
the condition for a family of vectors (ai,...,arf) to be a basis of M*^ is an open condition, be- 
cause it corresponds to det((ai, a^)) > where det is the determinant of the matrix. We 
apply the Gram-Schmidt process to the basis {h'{u),e2{z), ...,ed{z)) to obtain a direct orthonor- 
mal basis {h'{u)/\\h'{u)\\,V2{u, z), ...,Vd{u, z)) of for all u G Uz, such that the dependence 
of {h' {u)/\\h' {u)\\,V2{u, z), ...,Vd{u, z)) on u G Uz is of class C°° . We remark that the family 
{v2{u, z), ...,Vd{u, z)) is a direct orthonormal basis of Nr{h{u)) for all u G Uz. We have associ- 
ated with each z G r a neighbourhood Uz of Uz = h~^{z) in /, we can obviously suppose that 
Uz is an interval which is open in /. Since {Uz,z G r) is a covering of the compact /, we can 
extract a finite covering {Uj,j = l,...,n) from it. We can choose this family to be minimal, i.e., 
such that {Uj,j G {1, ...,n} \ jo) is not a covering of / for any jo G {1, ...,n}. We then reorder the 
(Uj,j = 1, ...,n) (keeping the same notation) by the increasing order of their left end point in I C M. 
Since the family (Uj) is minimal, each point of / belongs either to a unique set Uj, j G {1, ■■■,n}, or 
to exactly two sets Uj and f/j+i for j G {1, n — 1}. We denote by aj the middle of the non-empty 
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open interval Uj n f/j+i for j £ {1, ...,n — 1}, and by {h' {u)/\\h' {u)\\,V2{u, j), ...,Vd{u, j)) the direct 
orthonormal basis defined previously on Uj for j G {!,... ,n}. We want to construct a family of 
direct orthonormal basis {h' {u) /\\h' {u)\\, f2{u), fd{u)) oiW^ such that the function: 

: u e I ^ {h'{u)/\\h' {u)\\, f2{u), Uu)) 

is of class C°°. We have to define a concatenation of the {h' {u)/\\h' {u)\\,V2{u, j), ...,Vd{u, j)) over 
the different sets Uj. For u G [a, ai], we define 

i/j^u) = {h' {u)/\\h' {u)\\,V2{u,l), ...,Vd{u,l)) . 

Thus the function ip defined on [a,ai] is of class C°° . On [/inC/2 we have defined two different direct 
orthonormal basis {h' {u)/\\h' {u)\\,V2{u, j), ...,Vd{u, j)) for j = 1 and j = 2 that have the same first 
vector. Let (pi : Uir\U2 ^ 50d-i (M) be the function of class that associates to each u £ C/i n C/2 
the matrix of change of basis from {v2{u,2), ...,Vd{u,2)) to {v2{u,l), ...,Vdiu,l)). 

If bi is the right end point of Ui H U2, then 4>i is in particular defined on [ai, bi[. Let gi be a C°°- 
diffeomorphism from [ai, bi[ to [ai, oo[ which is strictly increasing (so 51 (ai) = oi) and such that all 
the derivatives of 51 at oi are null. Then (pi o g^^ is defined on [ai, +oo[ and all its derivatives at ai 
are equal to those of cfii. We then transform all the orthonormal basis {v2{u,j), ...,Vd{u, j)) of W^~^ 
for j > 2 and u > ai by the change of basis (pi o g^^ , and we denote the new direct orthonormal 
basis of R'^"^ obtained this way by {v2{u-,j), ■■■,Vd{u,j)). We then define tjj on ]ai,a2] by 

V'(n) = {h'{u)/\\h'{u)lv2{u,2),...,Vd{u,2)), 

and we remark that ip{u) still defines a direct orthonormal basis of M'^. The function xp is of class 
C°° on [a, 02], including at oi. We iterate this process with the family of basis 

{h' (u) /\\h' iu)\\,V2{u, j), ...,Vd{u, j)) , j = 2,...,n 

at a2, etc., finitely many times since we work with a finite covering of /. We obtain in the end a 
function 

i^oh-^ -.r^ SOd-i{R) 

which is of class C°°, and for all z £ r, the set of the points of that have for first coordinate in 
the basis ip o h~^{z) is exactly the hyperplane Nr{z). 

For each t = (t2, ■■■,td-i) £ {z £ R'^'^ \ d{z,0) < ??}, the set 

ft = {y £ \ 3z £ r , y has coordinates {0,t2, td-i) in the basis ip o h'^ (^z)} 

is a continuous path (even of class C°^) from a point in Xi to a point in X2, therefore 

nns nn ^ . 

Moreover, since d{S,r^ U F^) > 0, we obtain that 

nnsnn ^ 0. (lo) 

For each y £ tub(r, r/), there exists a unique Zy £ r such that y £ Nr{zy), so we can associate to y 
its coordinates {0,t2{y), ...,td{y)) in the basis ipoh~^{zy). We define the projection p of tub(r, r]) on 
^r{yi) that associates to each y in tub(r, r/) the point of coordinate {0,t2{y), ■■■,td{y)) in the basis 
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ipoh-'^{yi). Then p is of class as is if) oh ^. If z belongs to Nr{yi), and t{z) = (^2(2), ...,td{z)), 
then we know by equation (fTOl) that there exists a point on rj(^) that intersects S in fi. Moreover, 
rj(2) is exactly the set of the points y of tub(r, r]) whose image p{y) by this projection is the point 
z. Thus 

p{S n tub(r, 7]) net) = Nl!{yi) . 

Since tub(r, 77) is compact, p is a Lipschitz function on tub(r, ry), and so there exists a constant K, 
depending on p, hence on 17, r, ry, but not on 5, such that 

n'^-^{Sr\n) > W^-^S ntub{r,r])nn) > KW^-^ipiS ntuh{r,ri))) > Kaa-m'^^^ . 
This ends the proof of the positivity of (/>q when A(0) < 1 — Pc{d). 

Acknowledgment: The second author would like to warmly thank Immanuel Halupczok, Thierry 
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